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We derive an effective low-dimensional Hamiltonian for strongly interacting ultracold atoms in 
a transverse trapping potential near a wide Feshbach resonance. The Hamiltonian includes crucial 
information about transverse excitations in an effective model with renormalized interaction between 
atoms and composite dressed molecules. We fix all the parameters in the Hamiltonian for both one- 
and two-dimensional cases. 



I. INTRODUCTION 

Cold atomic gases strongly confined via optical lattice 
techniques along one or two dimensions and placed in an 
external magnetic field tuned near a Feshbach resonance 
provide the exciting possibility to study low-dimensional 
(low-D)strongly correlated physics in a controllable fash- 
ion [l|, . For weakly interacting atoms in a transverse 
trap, it is easy to get a low-dimensional Hamiltonian by 
projecting the atomic wave function to the lowest trap 
mode. However, the situation gets much more compli- 
cated for strongly interacting atoms. It is known that 
virtual excitation of the transverse trap modes during 
the atomic collisions could lead to interesting renormal- 
ization of the atomic scattering lengths in low dimen- 
sions 0, [IjQ- Furthermore, independent of the detuning 
from the Feshbach resonance, the atoms can always form 
dimers (two-body bound states) as their ground state 
configuration in a transverse trap [1, 0, H, H, 0, 13 ■ be- 
cause of the formation of dimers, the problem does not 
simply reduce to atomic scattering with a renormalized 
scattering length, but instead, we have to take into ac- 
count the real population of the atoms in transverse ex- 
cited levels, which remains significant even for an ex- 
tremely strong transverse trap [8|. 

In this paper, we construct the low-D (both ID and 
2D) effective Hamiltonian by taking into account the 
atomic population in the transverse excited modes. To 
describe the transverse excitations, one can introduce a 
dressed molecule Q , which is part of the dimer state ex- 
cluding the atomic population in the lowest transverse 
level. We observe that for realistic atomic densities, 
the structure of the dressed molecule is largely fixed by 
the two-body physics and thus almost density indepen- 
dent. This allows us to derive the renormalized interac- 
tion between the atoms and the dressed molecules, which 
gives the effective many-body Hamiltonian. This effec- 
tive Hamiltonian reproduces the correct two-body state, 
and from that, we also fix all the parameters in the Hamil- 
tonian at any detuning of the magnetic field. 



II. EXACT HAMILTONIAN 

We consider an atomic gas free in D dimensions {D = 
1, 2) and trapped by a (3— D)-dimensional harmonic po- 



tential of frequency to along the transverse direction. (A 
single well of the lattice potential can be well approxi- 
mated with a harmonic trap.) We use the conventional 
two-channel field theory for Feshbach resonance [lO|] , al- 
though one could also use a single-channel approach with 
an energy-dependent pseudopotential For atoms of 
mass m possessing internal states a ~ {ti i}j with con- 
tact interactions, the Hamiltonian is 

cr=T,i \ 4=1 / 

+ gb f (fr(^^\'i>l^ + h.c^ +Ub f d^r^-j^t^^^^ (i) 



where 5* (r) is the atomic field operator and $ (r) is the 
molecular field operator. Above, and throughout this pa- 
per, we use dimensionless quantities with all energies in 
units of hu! a nd all lengths in units of the trap length 
scale, at = y^h/muj. The dimensionless parameters are 
defined as follows: is the bare detuning (in units of 
huj), gb is the bare atom-molecule coupling rate (in units 

of hiuc^J"^), and J7{, is the bare background atomic scat- 
tering rate (in units of Huja^). The bare parameters are 
related to the physical ones via the standard renormal- 
ization relations: 
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where the subscript p denotes physical parameters, i^k = 
k^/2, and the integral is taken in three dimensions 
with an explicit energy cutoff Ec (which finally goes 
to infinity) imposed on two dimensions, so — 
V^E^/2'^/^7r. The physical parameters gp, Up, Vp are de- 
termined from the scattering data as Up — Anabg/at, 
gp = y^AnHcoWlabgl/athio, and i^p = ^co {B - Bq) /hio 
{pco is the difference in magnetic moments between the 
two channels), where we have assumed that the s-wave 
scattering length near resonance has the form = 
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, with Obg as the background scattering 
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length, W as the resonance width, and Bq as the reso- 
nance point. 

Expanding the field operators 'J (r) and $ (r) in terms 
of trap eigenmodes in the trapped dimensions and plane 
waves in the untrapped dimensions yields 
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where m indexes trap eigenmodes {rrii}, i — 1, . . . , 3 — 1?, 
and k denotes the wave vector in the untrapped dimen- 
sions {kj}, j = The operators flmko- and femk 
represent the corresponding atomic and Feshbach molec- 
ular modes, respectively. The plane wave energy Ck 
and the mode energy Sm are given by ek = ^i/-^' 

£m = X]i=i^ i^'^^' convenience we neglect the constant 
energy (3 — -D) /2 in as we will measure the two-body 
bound state energy with respect to the continuum thresh- 
old 3 — -D] . L is the dimensionless quantization length in 
the untrapped dimensions. The form factors appearing 
in Eq. ([3]) are given by 
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with 
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where Hn{x) is the Hermite polynomial. 

III. EFFECTIVE HAMILTONIAN 

The above Hamiltonian is extremely complicated and 
hard to solve directly. We expect it should reduce to some 
low-dimensional effective Hamiltonian when the trap con- 
finement is strong enough. In particular, we assume that 
the many-body energy scale, characterized by the den- 
sity dependent part of the chemical potential fip, is much 



less than the trap energy huj. This is equivalent to as- 
suming that the ID (or 2D) atomic density p <^ 1 mea- 
sured in units of a^^ (or a^"^), which is typically the 
case for realistic systems. Under this condition, when 
the atoms are far apart, they should stay in the trans- 
verse ground level to minimize the energy. The transverse 
excited levels get populated only when the atoms come 
close to strongly interact with each other. However, since 
p and there is no n-particle bound state with n > 3 
for two-component fermions, it is rare for three or more 
atoms to come close (two atoms can approach each other 
as there exists a bound dimer state at any detuning Vp 
from the Feshbach resonance [H). From this argument, 
we see that the atomic distribution in the transverse ex- 
cited modes is determined by the short-range physics, 
where the latter is fixed through the two-body solu- 
tion. The dimer state of the Hamiltonian H with bind- 
ing energy (relative to the continuum threshold, 3 — D) 
^E2B > and momentum q can be written in the form 
1^2) = where |0) denotes the vacuum and '^l^ = 



'^Oq'T^jnnk Bas -2ek-eq/2-e„ -£„ "'m,k-|-q/2,tXn, -k+q/2,i 

with the coefficients P and jymn given in Ref. Q , although 
the actual expressions are not important for the purposes 
of this discussion. We thus construct the dressed molec- 
ular modes djj with the same expression as ^2q' but 
excluding the m = n = term in the summation (cor- 
respondingly normalized). These modes c?^ capture the 
short-range physics and their structure should be basi- 
cally independent of the atomic density. The coupling 
between c?q and the atomic modes aoko- (with a simpli- 
fied notation as flko-) in the open channel can be approx- 
imated with a contact interaction since dq is tightly con- 
fined in space to a volume on the order of in the deep 
BCS regime and even smaller in the crossover and BEC 
regimes Q . The general effective Hamiltonian for c?q and 
flkcr then takes the form 
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where Afc (in units of hu)) is the relative detuning , on, (in 
units of fkua^^'^) is the coupling rate, and Vb (in units of 
hojaf) represents the remaining background interaction 
in the open channel. We introduce the physical parame- 
ters related to three bare parameters in H^g via a low-D 
renormalization analogous to Eq. 
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Note that the zero-point energy 3 — D appears expHcitly 
in the definition of V~^, otherwise there is an artificial 
infrared divergence. These definitions of the physical pa- 
rameters are justified as we will see below that they re- 
move exactly the ultraviolet divergence associated with 
the contact interaction. 

The effective Hamiltonian HcS should reproduce the 
same physics represented by the 3D Hamiltonian H when 
the system is near the ground state with the energy per 
particle close to (as /Zp ^ 1 in units of fuv). To de- 

termine the parameters in -ffefi, we match the exact two- 
body bound state obtained from the original H. Specifi- 
cally, we first determine the effective background scatter- 
ing Vp by matching the bound state physics in the deep 
BCS limit, where the population is entirely atoms in the 
lowest trap mode. Then, for general detuning, matching 
the binding energy and the bound state gives two con- 
straints which determine the remaining parameters, Af, 
and at,. Since the composition of the dressed molecule 
dq is a function of detuning Vp, or equivalently, of the 3D 
scattering length, so are these two parameters. 



IV. FIXING THE PARAMETERS 

All the two-body physics contained in a given Hamil- 
tonian are embodied in its T-matrix, defined by T{E) ~ 
Hi + HiG{E)Hi = Hi + HiGo{E)T{E), where Hi 
is the interaction part of the Hamiltonian H, G{E) = 
{E — H)~^ is the full two-body propagator, and Go{E) = 
{E — Hq)^^ is the free two-body propagator. Physically, 
the matrix element (2|r(i?)|l) is the sum of the di- 
rect process whereby a pair of atoms scatters from state 
|1) to state 1 2) and the indirect processes whereby the 
atoms scatter into an intermediate state, propagate with 
energy E, and then scatter into the final state. From 
its definition, it is clear that the T-matrix has simple 
poles where E is equal to a two-body bound state of the 
Hamiltonian. Furthermore, the residue of the diagonal 
matrix element {1\T (E) |1) at such a pole, E = En, is 

\{l\Hi\^n)\^ = \{l\H-Ho\^n)\^ = (En - e^f 

where ei is the noninteracting energy of state Thus, 
the residue of a given diagonal matrix element determines 
the fraction of the bound state in the given basis state. 
We shall use both of these facts in matching the two- 
body bound state properties of the exact and effective 
Hamiltonians. 

In the appendix, we derive the general two-body T- 
matrix associated with the Hamiltonian H in Eq. ([3]). 
For a pair of atoms asymptotically in the lowest mode of 
the trap, the corresponding diagonal T-matrix element 
(in units of aPfku) as a function of the two-body energy 
E (measured with respect to the continuum threshold, 
3 - L») is 

[TiE)]-'=jo^^[[UfiE)]-'-SpiE)) , (9) 
where 7000 = (27^)(^~'^/^ Uf{E) = Up ~ gl/{vp - E) 



and, from the normalization in Eq. (2), [E) = 

Uf {E) + U-^ = Ub- gll {vb -E) + U-\ and the func- 
tion 
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[13| . In the above we have used the gamma func- 
tion r (x) and the Hurwitz zeta function C,{s,x) = 

limAT^oo Y^ri=o ("^ + - + x)^"^ /(I - s). 

The two-body bound state corresponds to a pole of 
the T-matrix element above with [T (£'25)] ^ = at the 
binding energy |-B2b|, which gives the eigen-equation 



[Uf{E2B)] =Sp{E2B) 



(11) 



to determine E2b- Also, as discussed above, the atom 
pair population of the lowest trap mode is determined 
by the residue of the above T-matrix element at this 
pole. Since the bound state of the effective Hamilto- 
nian comprises only atom pairs in the lowest trap mode 
and dressed molecules, specifying the population in the 
lowest mode is sufficient, in conjunction with normaliza- 
tion, to determine the entire bound state of 7?off. Thus, 
to ensure that the effective and the exact Hamiltonians 
produce the same two-body bound state, we just need 
to obtain [T(i?)]~^ from both H and -ffeffj and require 
them to match to first order in _B — E2B ■ 

Following the same approach as in the Appendix, 
the corresponding T-matrix element from the effective 
Hamiltonian Hcs is obtained as 



(12) 



where [Vf (E)] ' = [Vp - a^/ {Xp - E) 
[Vb -all {Xb - E)] + V;"\ and the function 
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We require the effective background term alone to re- 
produce the two-body physics on the deep BCS side, 
where the dressed molecule population is negligible. 
Matching Eqs. ^ and (fT2|) in that region yields 



Vp-' - (2^) 



{3-D)/2 
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Cp) 
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where Gp = Hmi^^^^oo Sp (i?2s) - Cp {E2b)- In this way, 
we obtain an effective single-channel model on the deep 
BCS side which recovers the low-energy scattering mod- 
els of Refs. 0, |3:0|- However, the single-channel model 
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with a renormalized scattering length (oc T^) is not ad- 
equate near resonance or on the BEC side where the 
dressed molecule population becomes significant. 

To fix the parameters Xp and ap, we compare the T- 
matrix in Eqs. ([9]) and p2|) at general detuning, and 
require them to match for the zeroth and the first or- 
ders in expansion with {E — E2b) (which are responsible 
for reproducing the same binding energy and the bound 
state, respectively). After some algebra, we obtain 



A,] 



E. 



2B - 



<JpiE: 
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dE2 



Up^ ^ — {Sp — <7p)] 
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(15) 
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where 9£;,„ means d/dE\^ „ . Since the derivative of 
Eq. PU)) is not obvious, we write it below explicitly [12| |. 
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Thus, Eqs. (|14l[T5lfT6| along with the low-D renormal- 
ization procedure in Eq. ^ fix the parameters of the ef- 
fective Hamiltonian as functions of the two-body binding 
energy, E2B, which is connected to the physical detuning 
Vp through Eq. (jlip . The dctuning-dependent parame- 
ters are plotted in Fig. [1] across resonance, assuming the 
same typical 3D parameters as in Ref. The differ- 
ence in the sign of Up between ^Li and '^''K gives rise to 
quite different looking effective parameters. However, the 
relevant combination of parameters for mean-field calcu- 
lations, Vp^, is very nearly universal. We plot Vp^ (2/x) 
across resonance with fi = £^2b/2 -I- fJ-p for small values 
of /ip (/i has the meaning of the total chemical potential 
including the per-particle binding energy £^23 /2) . The 
slight difference between the ^Li and '"'K curves stems 
from the slight difference in binding energies as a func- 
tion of the 3D scattering length. In the low density limit 
lip 0, Vp^ can approach infinity in 2D on the BEC 
side where Ap — E2B, but in ID Vp^ is always attrac- 
tive since we have Ap > E2B at any detuning. The 2D 
resonance-like behavior {Vp^ {E) 00) around energy 
E ^ E2B should not be confused, however, with the con- 
finement induced resonance discussed in Ref. Q around 
energy E ^ 0. Due to the existence of the two-bod; 
bound state with E2B < 0, the resonance discussed in 
is not for the ground-state configuration of the system. 
Moreover, the relevant quantity for a typical many-body 
calculation is 1/Vp^ (E), which is continuous. 






^Li 2D 






---'"'K2D 






- - \i 1D 






"kid 











-a./a 



(a)Ap 



(b)a2 



t s 

2 





(c)2D 



(d)lD Vf 



FIG. 1: (Color online) (a)-(b): Detuning-dependent eflFec- 
tive parameters vs. inverse 3D scattering length, at /as — 
4TT/Uf{0). (c)-(d): Effective interaction (2/i) vs. at/a^. 
Solid lines are for ®Li and dashed lines are for ^"K. The curves 
correspond to /ip = 0, 0.05, 0.15, from left to right. 



other dimensions. Excited trap modes are important to 
the bound state physics, requiring the effective param- 
eters to assume a highly nontrivial magnetic field de- 
pendence. All the parameters in the Hamiltonian have 
been fixed from two-body considerations. This effective 
Hamiltonian can provide a starting point to understand 
the low-dimensional many-body physics when the system 
is near its ground state configuration (with the chemical 
potential close to the per-particle binding energy i?2s/2). 
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VII. APPENDIX 



V. SUMMARY 

In summary, we have derived an effective low- 
dimensional Hamiltonian for strongly interacting atomic 
gas trapped in one or two dimensions and free in the 



Below we derive the center-of-mass two-body T-matrix 
between atoms in arbitrary trap modes from the 3D 
Hamiltonian H in Eq. (1). We use |mnk) to denote an 
atomic pair state in trap modes m and n having relative 
momentum k, and |p) to denote a molecule in trap mode 
p. With the definition T = Hj + HiGqT, we evaluate 
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the propagator at energy E and obtain 



Substituting this form into Eq. ((20)) . we get 
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Substituting Eq. into Eq. ((THl), we get 
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where U^^ (E) is defined in the text. Note that we can 
separate out the final state dependence 

(mnk|T|m'n'k') = ^ Uf {E - Sp) 7mnpf^,„,^, (21) 



^m'n'k' ~ Tm'n'p + ^ {E £p') ^m'n'k' 
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The inner sum is equal to S {E — Sp) 5pip [1J|, where 
S [E) = Sp (E) ~ [/-I and 5*^ (E) is defined as in Eq. 
(fTO)l . so we can solve Eq. (|22]) to obtain 



jip 



and 
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For the special case of atoms asymptotically in the low- 
est mode of the trap, m = n = m' = n' = 0, we use the 
fact that 7oop Sp o to obtain Eq. ^ from Eq. ([24]) . 
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